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The structure of Lyapunov spectra for many particle systems with a random interaction between 
the particles is discussed. The dynamics of the tangent space is expressed as a master equation, 
which leads to a formula that connects the positive Lyapunov exponents and the time correlations 
of the particle interaction matrix. Applying this formula to one and two dimensional models we 
investigate the stepwise structure of the Lyapunov spectra, which appear in the region of small 
positive Lyapunov exponents. Long range interactions lead to a clear separation of the Lyapunov 
spectra into a part exhibiting stepwise structure and a part changing smoothly. The part of the 
Lyapunov spectrum containing the stepwise structure is clearly distinguished by a wave like structure 
in the eigenstates of the particle interaction matrix. The two dimensional model has the same step 
widths as found numerically in a deterministic chaotic system of many hard disks. 
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I. INTRODUCTION 



Chaos is characterized by a rapid expansion of a small 
initial error, and the Lyapunov exponent, which is de- 
fined as the time averaged exponential rate of the time 
evolution of infinitesimal perturbations of the dynamical 
variables, is introduced to express such characteristics 
quantitatively. The Lyapunov exponent is used to dis- 
cuss some statistical properties of many body systems 
(eg. the mixing property). It is also known that the 
Lyapunov exponents are connected to the amount of in- 
formation of the system (eg. the Kolmogorov-Sinai en- 
tropy), and to transport coefficients (eg. conductance 
and viscosity) . 

In general, the Lyapunov exponent depends on the di- 
rection of the infinitesimal perturbation of the dynami- 
cal variables at an initial time, so we obtain a Lyapunov 
exponent for each independent direction in the system. 
The sorted set of such Lyapunov exponents, is the so 
called Lyapunov spectrum, and has been the subject of 
study in many chaotic particle systems. For example, the 
thermodynamic limit of the Lyapunov spectrum, that is, 
that the spectrum retains its shape as the number of 
articles increases, is discussed using numerical evidence 
[j, random matrix approaches a periodic orbit 

approach JoJ and mathematical arguments |pqjLlJ| . Some 
works showed a linear behavior of the Lyapunov spec- 
trum in one dimensional models with nearest neighbor 
interactions although this conjecture is modified 

in weak chaotic systems fl2||l3f . The effect of the rota- 
tional degrees of freedom of molecules on the Lyapunov 
spectrum was investigated in a model consisting of di- 
atomic molecules, which showed an explicit separation of 
the rotational and translational degrees of freedom if the 
departure from sphericity is small enough |l4|,[l5). The 
Lyapunov spectra of systems in non-equilibriiim steady 
states were investigated by Refs. M n9[ with the dis- 
covery of the conjugate pairing rule for some systems 



with isokinetic thermostats |20|-p3|]. Avoided crossings 
and level repulsion in the Lyapunov spectrum, similar 
to the behavior of the energy levels in quantum chaotic 
system, was discussed in a mutually coupled map sys- 
tem |M . These works suggest the possibility of getting 
information about the system from the structure of its 
Lyapunov spectrum. The Lyapunov spectrum has also 
been of interest for the number dependence of the maxi- 
mum Lyapunov exponent in particle systems BJl9|, 25| , p6[ , 



and the Kaplan- York Lyapunov dimension 271- 30 [ 



One of the characteristics of the Lyapunov spec- 
trum, which has been shown recently in systems con- 
sist ing of many hard disks, is its stepwise structure 
0,^5 19.3l|j3^]. Such a stepwise structure appears in 
the distinct region of small positive Lyapunov exponents. 
Corresponding to this stepwise structure are the so called 
Lyapunov modes, a wave like structure in the tangent 
space of each eigenvector of a degenerate Lyapunov expo- 
nent, that is, for each stepwise structure. An explanation 
of the stepwise structure of the Lyapunov spectrum due 
to symmetries of the system was tried using two dimen- 
sional periodic orbit models ||. On the other hand, ex- 
planations for the Lyapunov modes were suggested using 
a random matrix approach for a one dimensional model 
p3| and more recently using a kinetic theoretical ap- 
proach p4| . The small Lyapunov exponents correspond 
to slow relaxation processes in the equilibrium state, so 
these results suggest the possibility of characterizing the 
macroscopic behavior in many particle systems from a 
microscopic point of view. 

This paper has two main purposes. First we consider 
the dynamics of the tangent space in dynamical systems 
with a random interaction between the particles. Dif- 
ferent from the other random matrix approaches for the 



Lyapunov spectrum ||||. 
els in time, we consider t 



|.33| which use discretized mod- 
;he continuous dynamics in time 
and derive a method to calculate the Lyapunov spectrum 
from a master equation technique. Using this method the 
Lyapunov spectrum is directly connected to the correla- 
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tions of the particle interaction matrix. As the second 
purpose we discuss the stepwise structure of the Lya- 
punov spectrum using this random matrix approach. We 
consider one and two dimensional models satisfying the 
total momentum conservation with periodic boundary 
condition. These models show a stepwise structure of 
the Lyapunov spectra only in a region of small positive 
Lyapunov exponents. We discuss the robustness of the 
stepwise structure against perturbations in the matrix 
elements of the particle interaction matrix. The long 
range interactions between the particles lead to a clear 
separation into a part of the spectrum exhibiting stepwise 
structure and a part changing smoothly with exponent 
number. We also find a wave like structure in the eigen- 
states of the particle interaction matrix in that part of 
Lyapunov spectrum containing stepwise structure. As a 
high dimensional effect we show that in the two dimen- 
sional model we get wider steps in the Lyapunov spec- 
trum, such as steps consisting of 4 and 8 degenerate ex- 
ponents, rather than simply 2 exponents as in the one 
dimensional models. The step widths of 4 and 8 expo- 
nents are actually observed in the numerical simulation 
of the many hard-core disk system [pTJ . 



II. RANDOM MATRIX DYNAMICS OF THE 
TANGENT VECTOR 



We consider a Hamiltonian system whose state at time 
t is described by the 2A^ dimensional phase space vec- 
tor T(t) = (q(t),p(t)) T with the coordinate vector q(£), 
the momentum vector p(t) and the transpose operation 
T. The dynamics of the phase space is described by the 
Hamilton equation 



dT(t) _ T dH(T(t),t) 



dt 



= J 



dT(t) 



(1) 



with the Hamiltonian H(T(t),t). Here we allow an ex- 
plicit time-dependence of the Hamiltonian, and J is the 
2N x 27V matrix defined by 



J : 



/ 
-I 



(2) 



where / is the N x N identical matrix and is the N x N 
null matrix. It should be noted that the matrix J sat- 
isfies the conditions J T =_—J and J 2 = —I with the 
2N x 2N identical matrix I. 

The dynamics of the tangent vector ST(t) as a small 
deviation of the phase space vector is described by the 
linearized equation 



dST(t) 
dt 



JL(t) ST(t) 



(3) 



of the Hamiltonian equation (|l]) where L(t) is defined by 
_ d 2 H(T(t),t) 



L(t) 



dT(t)dT(t) 



(4) 



The matrix L(t) is symmetric: L(t) — L(t) T . At this 
point we note that an equation equivalent to Eq. (j^) 
can also be obtained for system whose dynamics is not 
derivable from a Hamiltonian. 

Equation (|J) is formally solved as 



ST(t) = M(t, t )5T(t ) 
where the matrix M(t,to) is defined by 



M(t,t ) =T exp 



J / ds L(s) 

to 



(•5) 



(G) 



with the positive time-ordering operator T with the lat- 
est time to the left. Equation (H) leads to 



M^to)- 1 = T exp 



J / ds L{s) 



= -JM(t,t ) T J, 



(7) 



with the negative time-ordering operator T with the lat- 
est time to the right, so that the Hamiltonian phase vol- 
ume is preserved, namely |Det{Af (i, to)}\ = 1. 

For simplicity, we consider the case of zero magnetic 
field, so the Hamiltonian H(T(t),t) is represented as 
H(T(t),t) = |p(i)| 2 /(2m) + V(q(t),t), where m is the 
mass of particle and V(q(t), t) is the potential energy at 
time t. In this case the matrix L(t) is given by 



L(t) = 



-R(t) 
I/m 



(8) 



where R(t) = (i? a( g(t)) is the N x N symmetric ma- 
trix defined by -d 2 V{q(t),t)/dq{t)dq{t). The effects of 
particle interactions and fixed scatterers are taken into 
account through the matrix R(t). 

Now we consider the case where the particle inter- 
actions or collisions with fixed scatterers occur ran- 
domly enough, so that the matrix elements R a p(t), a = 
1, 2, • • ■ , N, f3 = 1, 2, • • • , N can be regarded as Gaussian 
white randomness in the sense of 

(Rai/3i(tl)Ra 2 2 (h) ■ ■ ■ -Ra 2 n-i/32n-i(*2n-l)) =0, (9) 
(-R«i/3i (ti)R a2 f3 2 (£2) ' ' ' Ra 2n 2n (*2n)) 

^ ^ ^ a 31 031 a 32 ^ a 33 @J3 a J4 
Pd 

~^ a 32n-l ^J2n-1 Q J2n ^32n 

xS(t h - t h )5(t h - t u ) ■ ■ ■ 5{t 32n _ t ~ t j2n ) (10) 

for any integer n, where we take the sum over only the 
permutation P d : (1,2,- •• ,2n) -> (J1J2, • ■ • Jin), and 
the bracket (• • •) means the ensemble average over ran- 
dom processes. Here Djkin is a 4-th rank tensor and is 
assumed to be constant in the rest of this paper. The 
tensor Djkin satisfies the condition 



Injk 



D 



jkln 1 



(11) 



2 



because the equation Djki n S(s — t) = (Rjk(s)Ri n (t)} 
= (Ri n (t)Rjk(s)) = Di n jk5(s — t) is satisfied at arbitrary 
times s and t. The symmetric property of the matrix 
R{t) also imposes the conditions 



D 



jknl 



— Dkjln — D 



jkln 



(12) 



for the tensor Dj^in- 

Under the Gaussian white conditions (|J) and (|l0) for 
the matrix R(t), the dynamics described by Eq. ([3) for 
the tangent vector can be regarded as a stochastic pro- 
cess. Now we consider the description of this stochastic 
process by a master equation for the probability density 
p(ST, t) for the tangent vector space in which p(ST, t)dST 
is introduced as the probability of finding a tangent vec- 
tor ST in the region (ST, ST + dST). By applying the 
Kramers-Moyal expansion to the dynamics (|3|) with the 
randomness given by Eqs. (|^) and ( [To| ) we obtain 

dp(5T,t) 



dt 



N 

E 

a = l 

N N 



6 Pa dp{ST,t) 
m dSq a 

JV JV 



E E E E Q5 d5pv 



a=l ;j=l )l = l V=l 



(13) 



Equation (U3) is a master equation for the tangent vec- 
tor ST = (7q, 5p) T with <5q = (Sqi, 5q2, ■ ■ ■ , SqN) T and 
Sp = (Spi, Sp 2 , • ■ ■ , Spn) t , especially the type of the 
equation called Fokker-Plank equation because its right- 
hand side includes up to the second derivative of the 
probability density with respect to the variables. The 
derivation of Eq. (13) is given in Appendix |A|. We should 
notice that the first term of the right-hand side of Eq. 
(H) has the same form as the corresponding term in the 
Liouville equation describing non-interacting particle sys- 
tems, because the dynamics of the phase space point and 
the tangent space point coincide in the non-interacting 
particle system. The characteristics of the system are 
introduced through the 4-th rank tensor D a ^0 v and a 
boundary condition for a solution of Eq. dl3). 



III. LYAPUNOV SPECTRUM IN THE RANDOM 
MATRIX APPROACH 

In this section, using the Fokker-Plank equation Ji~3| ) 
for the tangent vector we investigate the absolute values 
of the tangent vectors, whose asymptotic behaviors lead 
to the Lyapunov exponents. 

(0/ 



We introduce the quantities T) fe (i), I = 1, 2, 3 as 



r$(t) = (6 qj 6q k ) t 
Tf k \t)^((Sq 3 Sp k + Sq k Sp 3 ))j2 
rf k \t)^(5p 3 5 Pk } t 



(14) 
(15) 
(16) 



where the bracket (• • -)t with the subscript t means the 
average over the probability density p(ST, t) at time t: 



(• ■ -)t = J dSTp(ST, t) ■ ■ ■. These quantities consist of el- 
ements of the averaged matrix (<5r<5r T ) t , and satisfies 
the condition 

T W(t) = Tg(t). (17) 
Equation ( |l3| ) connects these quantities as 



dt 

*(2) 



dt 

drf k \t) 

dt 



m 

JV JV 



= EE D ^ T $w 



(18) 
(19) 
(20) 



where we assumed the probability density p(ST,t) to be 
zero at the boundary of the tange nt sp ace, and used Eqs. 
(|ll]), © and (O) to derive Eq. @. 

Now we introduce a 4-th rank tensor Tj k i n satisfying 
the condition 



JV JV 

E E TajkpTpnla — SjlSkr, 
a=l /3—1 



(21) 



so that we obtain 



JV JV JV N 



E E E E Tj a pkD afiV pT nI/f/t i — AjkdjiSkn (22) 

a=l (3=1 ii=l v-1 

with the real 2nd rank tensor Ajk satisfying the condi- 
tion A j3 > @. (Concerning Eqs. @ and (§§), for 
example, the existence of the real tensor Ajf. satisfying 
Eq. (22) simply comes from the fact that by the condi- 
tion (12) we can regard the quantity Dj^in as the matrix 
element 2->7i(j,n)7 2 (fc,0 of the N 2 x ./V 2 real symmetric ma- 
trix V = (^ , 7l (j, n ) 72 (fe,;)) where j n = 7„(j, k), n = 1, 2 are 
functions from j G {1, 2, • • • , N} and fee {1, 2, • • • , N} to 
7„ G {1, 2, • • • , N 2 }.) We will discuss an example of the 
tensors Tjkin and Ajk in the next section. By using the 

tensor TjMn we transform the quantities Tv (t) to xS(i) 



JV JV 
c*=l /3=1 



(23) 



The inverse transformation to derive the quantity (t) 
from the quantity T„«(i) is simply given by 



JV JV 
a=l (3=1 



(24) 



noting the relation @. Especially we should notice 
the relations (|£q| 2 ) t = ELi f ««(<) and (I^P| 2 >t = 
Ea=i ^ai(t) under the condition J2a=i T jaak = <5jfe- 
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We introduce the positive (or zero) Lyapunov expo- 



nents A,- as 



Xj = lim — In 



(25) 



namely as the exponential rate of the time evolution of 
the average of magnitude of the infinitesimal deviation of 
the phase space orbit in the long time limit. As shown in 
Appendix [l|, these Lyapunov exponents are simply given 
by 



A, 



(2m) 



1/3 



(26) 



using the quantity Ajj introduced in Eq. (j22|). Equa- 
tion (p6| ) is the key result of this paper. This equation 
connects the Lyapunov exponents directly with the ten- 
sor Djkin representing the strength of correlation of the 
particle interactions and potentials given by fixed scat- 
terers, and also shows the fact that in the system de- 
scribed by the random matrix dynamics the Lyapunov 
exponents are independent of the initial condition like 
in the deterministic chaos. It should be noted that 
the Lyapunov exponents Xj given by Eq. (|2^) are the 
same with the quantities derived from the equations 
lim t ^ oo (2t)- 1 ln[fg ) (t)/fg ) (0)],/c = 2,3 in this ap- 
proach (See Appendix |b|). 

We sort the Lyapunov exponents Ai, A2, Ajv so 
that they contain an decreasing (or equal) sequence, and 
introduce the set {A' 1 ', A^, • • • , A^} satisfying the con- 
dition 

{A[ 1 ],A[ 2 ],---,A[ Ar l} = {A 1 ,A 2 ,---,A J v} (27) 



AW > AM >.-->A^. (28) 

This sorted set of the Lyapunov exponents is called the 
Lyapunov spectrum. 



IV. A SIMPLIFICATION 

Before considering Lyapunov spectra in concrete mod- 
els using the formula given in the previous section, we 
discuss an assumption to simplify model calculations. 
We consider the case where all matrix elements Rj k (t), 
j = 1, 2, • • • , TV and k = 1, 2, • • • , N have the same time 
dependence, namely 



R jk (t)=r(t)A jk 



(29) 



where r(t) is a normalized Gaussian white randomness in 
the sense of (r(ii))r(i 2 ) ■ • ■rfon-i)) = and (r(ti)r(t2) 
■■■r(t 2n )) = Ep d H t n- t j2) s ( t j 3 - t u) ■•• ^(*j2 B -i-*32„) 
for any integer n and A = {Ajk) is a time-independent 
N x N matrix. The matrix A must be a symmetric ma- 
trix: 



.4 



A 



(30) 



because the matrix R(t) is symmetric. 

The assumption (E9t) simplifies our considerations, be- 
cause in this case the tensor Dj k i n is represented as the 
multiplication of the matrix element Aj k with A\ n . Thus 



D 



jkln 



A ik A 



jk A ln 



(31) 



and the conditions (Jllj) and (|12|) for the tensor Dj k i n 
are automatically satisfied under the condition (]30|). The 
condition ( ^0| ) also implies that the matrix A is diagonal- 
izable using an orthogonal matrix U = (Uj k ) satisfying 
U T U = UU T = I, namely (U T AU) jk = ajS jk with real 
eigenvalues a;, / — 1, 2, • • • , N of the matrix A. The 4-th 
rank tensor Tj k i n is constructed as 



T. 



jkln 



UkjUl r , 



(32) 



which satisfies the conditions (pi"]), (53) and 



53a=l Tjack — 3jk- 

represented as 



Here the 2-nd rank tensor Aj k is 



A jk = aja k (33) 

which satisfies the conditions Ajj > 0. Equations ( pjjj ) 
and d33| ) lead to the expression of the Lyapunov expo- 
nents as 



A; 



2m 



2/3 



(34) 



After all under the assumption (|29j) the calculation of 
the Lyapunov spectrum is attributed to the eigenvalue 
problem of the matrix A. 



V. ONE DIMENSIONAL MODELS AND STEPS 
OF THEIR LYAPUNOV SPECTRA 



In this section we consider simple one dimensional 
models and calculate their Lyapunov spectra using the 
formula given in Section III under the assumptions dis- 
cussed in Section IV. We are generally interested in 
which ingredients of the model system lead to particu- 
lar features of the Lyapunov spectra. We are especially 
interested in models which satisfy total momentum con- 
servation and show stepwise structures in their Lyapunov 
spectra. Numerically the observation of Lyapunov modes 
is associated with systems with a stepwise structure of 
their Lyapunov spectrum, such as the many hard disk 
system of Refs. ||||^,||^|. 

We construct a model consisting of N particles in a one 
dimensional space. In this case the off-diagonal matrix 
element Aj k , j ^ k represents the strength of the interac- 
tion between the j'-th particle and the fc-th particle. The 
diagonal matrix element Ajj is determined by the total 
momentum conservation which imposes the condition 



A' 



5>» 

fc=i 



(35) 
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for the matrix A. Equation ( |35[ ) is derived from Eq. (|2 

and the equation yi^l-, Rjk(t) = satisfied at any time 
t. Using the formula (pj) this condition implies that 
the iV-dimensional vector, whose components are equal, 
is an eigenstate of the matrix A corresponding to the 
eigenvalue 0, so we obtain 



\W = 0, 



(36) 



namely there is a zero Lyapunov exponent corresponding 
to the conservation of total momentum. 



where Xji 3 = 1> 2, • • • , N + 1 are real constants satisfy- 
ing the condition xn+i — Xi- I n this case we can calcu- 
late the eigenvalues of this matrix A to first order in the 
small parameter e, and obtain the Lyapunov exponents 

A„ = An + eXn* 1 + 0(e 2 ) in the expanded form by the 



parameter e with the quantities 



2A 



(0) 



3w 



(Xo + |Xil ) 



2A 



(0) 



,(1) _ *"N-j 



3u) 



(Xo - |Xil ) 



(41) 
(42) 



A. One dimensional model with a stepwise structure 
of its Lyapunov spectrum 

As a first step, we consider the case where each parti- 
cle interacts only with its nearest neighbor particles with 
the same strength of interaction. We impose periodic 
boundary conditions, namely that the particles are on 
a one dimensional ring structure. This situation is de- 
scribed by the matrix A — A^ = (A 1 '?) defined by 



1(0) 



-26 



jk 

+5 



7 U + l)k 



j(k-N+l) 



i) + 



U-N+l)] 



(37) 



with a (non-zero) real constant cu. We can calculate the 
eigenvalues of the matrix A^ (See Appendix |c| for the 
calculation details.), and obtain the Lyapunov exponents 

A,, 



x (o) 



as 



|w| ( 2irn\ 
— 1 - cos — — 
m V N I 



2/3 



It is important to note that the Lyapunov exponents 
given by Eq. ( ^8| ) satisfy the conditions A^ = and 



x (o) _ AO) 



(39) 



for j < N/2, namely the Lyapunov spectrum has de- 
generacies. In other words the spectrum has a stepwise 
structure. The maximum Lyapunov exponent is given by 
(2|w|/m) 2 / 3 when N is even or {|w|[l + cos(7r/A)]/m} 2 / 3 
when N is odd. 



B. Robustness of the stepwise structure of the 
Lyapunov spectrum to a perturbation 



Next we consider the case where each particle interacts 
with its nearest neighbor particles by slightly different 
interaction strengths. This situation is described by the 
matrix A = A^ + eA^ with e a small parameter and 
the matrix A^ defined 



A fk = ~(Xj + Xj+i)$jk + XjSj{k+i) + XkSk(j+i) 



-Xi [Sj(k- 



N+l) 



J k (j-N+1)\ 



(40) 



in the case of j < N/2, where Xj 1S defined by 



Xi 



1 / AiriajX 

a=l v ' 



(43) 



(See Appendix ^| for the derivations of Eqs. ( pd] ) and 



) 



Here we numbered so that we obtain w(A^ — 

> 0. It may be noted that the quantity Xo i n 
flU ) and ([l2|) is simply the arithmetic average of 
the coefficients Xj, j = L 2, • • • , N: xo = N^ 1 J2f=i Xj- 
2A]^Lxo/(3w) in the case of the 



Eqs. 



We also obtain X^, 2 ~ 
number N to be even, and A 



(i) 

N 



for the first order 
correction to the non-degenerate Lyapunov exponents. 

Equations (|4l]) and (|4|) tell us about robustness of 
the stepwise structure of the Lyapunov spectrum appear- 
ing in the non-perturbed system against the perturba- 
tion. To discuss this point we consider the deviation 
(38) |Aj — Ajv-j | of the two points in the j-th step of the Lya- 



punov spectrum, which is given by 



|A, - A 



N-j\ 



£ Xj 



X, +0(e 2 ) 



(44) 



in the case of j < N/2. Equation ( f44| ) implies that the 
degeneracies, namely the stepwise structure of the Lya- 
punov spectrum is removed by the perturbation sA^ 1 ' , 
but their deviation are small in a region of small Lya- 
punov exponents rather than in a region of large Lya- 
punov exponents, as far as the quantity Xj is almost j- 
independent. In other words, this consideration suggests 
that the stepwise structure of the Lyapunov spectrum is 
robust in the region of small Lyapunov exponents. It 
should be emphasized that this is consistent with numer- 
ical results of the Lyapunov spectrum in a many hard 
disk model, which shows the stepwise structure of the 
Lyapunov spectrum only in a region of small Lyapunov 
exponents . By using Eqs. @ and @ 

we can also discuss a perturbational effect in a global 
shape of the Lyapunov spectrum. Eqs. © and @ 
lead to the shift of the j-th step of the Lyapunov spec- 
trum as 



Aj + X N -j 



A 



(o) 



2xo 
3cj 



A, 



■0(e 2 



(45) 



5 



in the case j < N/2. Equation ( f45j ) implies that in the 
case of Xos/ w > (xqe/uj < 0) the perturbation eA^> 
makes the slope of the Lyapunov spectrum more (less) 
steep than in the non-perturbed case. 



^ ] = E[-(^ ] +^)W ] 



j-°j(k+i) +v k 

(46) 



i [/] c , [I] X 

+ a j-N+l°j(k+N-l) + a k _ N+l k (J+N-l) 
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FIG. 1. Lyapunov spectra of a one dimensional system with 
nearest neighbor interactions in the case of randomly chosen 
matrix elements. The inset: Enlarged graph of a small Lya- 
punov exponent region. 

Now we discuss the robustness of the stepwise struc- 
ture of the Lyapunov spectrum against perturbation in a 
little different way, in which the fluctuations of the ma- 
trix elements of A are large enough so that the above 
first order perturbed discussion is no longer correct. We 
consider a one dimensional system described by the ma- 
trix A = A^ given by Eq. (Efl), in which the quantities 
Xj, j = 1,2, •••,N are chosen randomly from the re- 
gion (a, b). Figure Q is a Lyapunov spectrum normalized 
by the maximum Lyapunov exponent in such a system 
consisting of 100 particles (N — 100). Here, we chose 
the region (a, b) as a = uj/2 and b = 3u>/2, and took 
the arithmetic average over this randomness for the nor- 
malized Lyapunov spectrum. Although the magnitude 
b — a (= u) of fluctuations in the matrix elements of the 
matrix A is the same scale as the averaged magnitude 
(a + b)/2 (= u>) of the matrix elements of the matrix A, 
we can still recognize some steps in the Lyapunov spec- 
trum in the region of small Lyapunov exponents. It is 
interesting to note that in such a case the global shape of 
the Lyapunov spectrum is rather close to a straight line, 
like discussed in Refs. j| ||. 



C. Effect of long range interactions and wave like 
structures in eigenstates 



In this subsection we consider the effects of long range 
interactions between particles on the Lyapunov spec- 
trum. The effect of the long range interactions between 
particles can be taken into account using the matrix 
A = A^ = {A [ ™ l] ) defined by 



with real constants ay , where ni(< N/2) is the length of 
the interactions and its value means that particles can in- 
teract with up to their ni-th nearest neighbor particles. 
As shown in Appendix O, we can obtain an analytical 
expression for the Lyapunov spectrum derived from the 
matrix (|4^) if all of the non-zero matrix elements erj'' , 
j = 1, 2, • • • , N, I = 1, 2, • • • , rii take the same value. 




j 

FIG. 2. Lyapunov spectra of one dimensional systems with 
different length n; of long range interactions in the case of 
N — 100. The three graphs correspond to the cases n; = 10 
(the circle dots), n; = 15 (the triangie dots) and n; = 20 (the 
square dots), respectively. The inset: Enlarged graphs of a 
small Lyapunov exponent region. 

In this subsection we consider the case where the ma- 
trix element er^ is chosen randomly, like in the model 
discussed in the end of the previous subsection. Figure |2| 
is the Lyapunov spectrum normalized by the maximum 
Lyapunov exponent in the system described by the ma- 
trix A = A^ with the quantities chosen randomly 
from the region (uj/2, 3a>/2) for the case of ni = 10, 15 
and 20 in the system consisting of 100 particles (N=100). 

(This randomness in the quantities oy^ is adopted to draw 
all of the figures hereafter in this subsection.) In these 
graphs we took the arithmetic average of the Lyapunov 
spectra over the randomness of the quantities . This 
figure shows that the long range interactions separates 
the Lyapunov spectrum clearly into a part exhibiting 
stepwise structure and a part changing smoothly. (The 

randomness of the quantities cn is not essential for this 
separation in the Lyapunov spectrum. It mainly plays 
the role of smoothing the Lyapunov spectrum in the re- 
gion of large Lyapunov exponents.) The stepwise struc- 
ture of the Lyapunov spectrum appears only in a region 
of small Lyapunov exponents, as suggested in the previ- 
ous subsection. 
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FIG. 3. Dependence of the long range interaction length n; 
on the Lyapunov exponents A' J ' in the case of N = 100. The 
graphs with the circle dots, the triangle dots, the square dots, 
the open circle dots and the open triangle dots correspond to 
the case of j = 99, 97, 95, 93 and 91, respectively. 



The stepwise region of the Lyapunov spectrum is 
smaller and smaller as the interaction length ni of par- 
ticles is longer and longer. (As a limit of the long range 
interactions, as shown in Appendix we can easily show 
that in the system, where each particle interacts with all 
the other particles with the same strength, all of the pos- 
itive Lyapunov exponents take the same value.) Besides, 
the heights of the steps of the Lyapunov spectrum are 
approximately proportional to the length ni of the inter- 
actions. These characteristics are clearly evident in Fig. 
|3|, showing the length ni dependence of the j-th posi- 
tive Lyapunov exponents A^'l , j = 91, 93, 95, 97, 99 in the 
case of 100 particles. In these graphs we also took the 
arithmetic average of the Lyapunov exponents over the 

randomness of the quantities . This figure shows that 
the values of small Lyapunov exponents increase with the 
interaction length ni as far as they are in the region of 
the steps of the Lyapunov spectrum, and if the interac- 
tion length ni is bigger than a critical value then their 
ri; dependences are slowed down, meaning that they are 
in the region of the Lyapunov spectrum that is changing 
smoothly. 

It may be interesting to investigate the difference be- 
tween the two parts of the Lyapunov spectrum from the 
point of view of the eigenstate. Figure [| is the real eigen- 
state of the matrix in the case of ni — 15 and 

N = 100 with randomly chosen quantities Oj\ (Note 
that we did not take the average over the randomness 
of the quantities to draw this figure.) The graphs 
(a), (b), (c) and (d) are the eigenstates corresponding to 
the Lyapunov exponents for j = 99, 98, j — 97, 96, 
j — 93,92 and j = 91,90, respectively. These Lyapunov 
exponents correspond to the triangle dots surrounded by 
the broken lines in the inset of Fig. ||. We can see clear 
wave like structures of approximately sinusoidal type in 
the graphs (a) and (b), which correspond to the Lya- 
punov exponents composing the steps in the Lyapunov 
spectrum. The wavelength of the waves in the graph (b) 
is half of the wavelength of the waves in the graph (a). 



On the other hand, we cannot recognize such a wave like 
structure in the graphs (c) and (d), which belongs to the 
Lyapunov exponents in the part of the Lyapunov spec- 
trum which is changing smoothly. 
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FIG. 4. Eigenstates of matrix A^" 1 ^ in the case of ni = 15 
and N = 100: the graphs (a) (j = 98,99), the graph (b) 
(j = 96,97), the graph (c) (j = 92,93) and the graph (d) 
(j = 90, 91) for the Lyapunov exponent A' J '. These Lyapunov 
exponents correspond to the triangle dots surrounded by the 
broken lines in the inset of Fig. 



It should be noted that the wave like structure of the 
eigenstates already appears even in the system described 
by the matrix A = A^ defined by Eq. ( p7\ ) with nearest 
neighbor interactions of a constant strength. Therefore 
the point is that such a wave like structure of the eigen- 
states is not destroyed by the long range interactions and 
the random interaction strengths only in the eigenstates 
which correspond to the small positive Lyapunov expo- 
nents in the steps of the Lyapunov spectrum. 

One may regard the long range interactions discussed 
in this subsection as a kind of high dimensional effect, 
because if particles move in a two or three dimensional 
space then they can interact with more than two parti- 
cles, even in hard-core interactions. However, some high 
dimensional effects, for example, the total momentum 
conservation in each orthogonal direction and so on, are 
still missing in the discussions in this subsection. In the 
next section we consider high dimensional effects more 
explicitly. 



VI. TWO DIMENSIONAL MODEL WITH A 
TRIANGLE LATTICE STRUCTURE 



In the one dimensional models discussed in the pre- 
vious section, the steps of the Lyapunov spectrum are 
caused by the periodic boundary conditions, and all of 
them consists of 2 only points. In this section we consider 
a two dimensional model, and show that it is possible to 
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get wider steps than in the one dimensional model as a 
high dimensional effect. Specifically, we concentrate on a 
model which has steps consisting of 4 and 8 points, which 
are actually found in a square system consisting of many 
hard disks 

In the two dimensional system, the position of each 
particle is specified by a two dimensional vector q^, 

j = 1, 2, • • • , N, where N is the number of the parti- 
cles and N = 2N. In such a case the matrix A is a 
(2N) x (2N) matrix, and can be represented as a block 
matrix consisting of 2 x 2 matrices B^ k ' , j = 1, 2, • • • , N 
and k = 1, 2, • • • , N. Here B^ = (B<f$) is given by 



A 



(2i-l)(2fe-l) A(2j-l)(2k) 



A 



(2j)(2fc-l) 



A 



(2j)(2fe) 



(47) 



corresponding to the matrix — d 2 V/dqjdq k , and its com- 
ponents represent the strengths of the interactions be- 
tween components of positions of the j-th particle and 
the A:-th particle. The matrix B^ k ^ is symmetric and 
the matrix A is also symmetric in the sense of the block 
matrix, namely 



B 



j'k' 



= B 



W 

k'j' 



B 



j'k' ■ 



(48) 



The diagonal block bCm') is determined by the condition 
of the total momentum conservation 



N 



fe=l 



which comes from Eq. 



(49) 



and the relation 



SfeLi ^V/dlj^Qk = an( i i s simply the two dimen- 
sional version of the condition ( |35| ) . This condition J49] ) 
implies that the Lyapunov exponents include at least two 
zero components: 



X [N-1] = X [N] = 



(50) 



corresponding to the conservation of the total momen- 
tum. 

To construct a two dimensional model it is convenient 
to use a lattice picture, because in the approach of this 
paper, under the assumption discussed in Section [V, the 
model is specified by pairs of interacting particles and 
their interaction strengths. In such a lattice, each lattice 
site corresponds to a particle and a connection between 
sites means that particles on those two sites can inter- 
act with each other as nearest neighbor particles. Now 
we consider an equilateral triangular lattice from such a 
point of view. This situation is motivated by the fact that 
in Ref . |3|] the system consists of many disks which pack 
in a two dimensional space as a hexagonal close-packed 
structure in the high density limit. The triangular lattice 
has three directions to connect the sites, and arrange the 
triangle lattice so that it fits within a square where one 
of sides of the square is parallel to one of the three direc- 
tions of the triangle lattice, as shown in Fig |B|(a) where 



such parallel lines are horizontal. For simplicity we as- 
sume that the number of lattice sites in each horizontal 
line is equal and is given by iVi(> 1). We put the num- 
ber of such lines as N 2) which should be roughly given 
by JV 2 « (2/V3)N v 



(b) 
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FIG. 5. (a) The triangular lattice system with a square 
shape. JVi is the number of the particles in a horizontal line, 
and N2 is the number of the horizontal lines, (b) Number- 
ing of the particles and the boundary condition in the part 
surrounded by the square of broken line in the figure (a). 



Next, we introduce the boundary condition for this tri- 
angle lattice model. For this purpose we assign numbers 
1, 2, • • • , JVi to the lattice sites (namely particles) in the 
first horizontal row of particles from left to right, and 
numbers N\ + 1, Ni + 2, • • • , 2Ni in the second row and 
so on until the last row numbered (N2 — l)Ni + 1, (A^2 — 
l)7Vi + 2, • • • , N 2 Ni, as shown in Fig. |(b). We define 
the number /23-1 as the number of the particle which do 
nearest neighbor interactions with the [2(j — l)Ni + 1]- 
th particle, and the number f 2 j as the number of the 
particle which do nearest-neighbor interacts with the 
[2(j - l)Ni + l]-th, [(2j - l)JVi + l]-th and [2jJVi + l]-th 
particles (j — 1, 2, • • • , Int{A^2/2}). Here, Int{x} means 
to take the integer part of x for any real number x. We 

also define the numbers and ^ as the number of 
particles which do nearest neighbor interactions with the 
j-th particle through the upper-left line and the upper- 
right line, respectively (j = 1, 2, • • • , Ni). (See Fig. S(b) 

for the definitions of the numbers fj and .) The num- 
ber fj specifies the interaction between the left side and 
the right side of the square, and the numbers g^ spec- 
ifies the interaction between the upper side line and the 
lower side of the square. In this section we consider the 
case that these numbers are given by 



fj = jNi 



5l (±) = (Ni - 1)N 3 + min 
where is defined by 



(51) 



l jk 



> 



} (52) 
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hff=j±Int{~[N 2 - 



1 ± 1 



(53) 



The boundary conditions ( |5l| ) and fl52|) imply that every 
line of the three directions in the triangle lattice is peri- 
odic. We assume the condition ^= g\~' so that every 
particle has the six nearest-neighbor particles. 

Before giving the matrix A of the triangle lattice model 
including a long range interaction, we consider the ma- 
trix A — Ao including only the effect of the nearest 
neighbor interactions of a constant strength. We de- 
fine the 2Ni x 2iVi matrices Ci, I = 1,2,3,4, each of 

(ik) 

which is a block matrix consisting of 2 x 2 matrices Cf , 
j = l,2,---,N 1 ,k = l,2,---,N 1 ,by 



C 



(Jk) 



°j(k-N 1 +l) 



h(j-N 1+ i)] (54) 



Cj* = B [Sjk + $j(k+l) + $j(k-N 1 +l)] 



C W=B[5 jk + 5 k 



fc(i+i) 



3fc(j-JVi + l)J 



(55) 



(56) 



C i k) = B { S k [g ^-(N 1 -l)N 2 ] + 6 k[ g ' j ^-(N 1 - 



where B is a 2 x 2 matrix. The matrix Aq is introduced 
as the block matrix defined by 



A = 



( °1 


^2 








• 


• 


c 4 \ 


cj 


Ci 


c 3 





• 


• 








ci 


Ci 


c 2 















ci 


Ci 


c 3 















ci 




















Cx 


Cp 













ci 


CiJ 



(58) 



with the 2Ni x 2Nx null matrix 0, where v is 2 (3) if N 2 
is an odd (even) number. 

The matrix A = A^'l of the triangle lattice model in- 
cluding the effect of the long range interactions up to the 
n;-th nearest neighbor interactions is simply given as fol- 
lows. We consider the matrix A' given by the matrix Aq 
except that the matrix B in the matrix A are replaced 
by the 2x2 matrices whose matrix elements are positive 
constants. By using such a (2N1N2) x (2N1N2) matrix 
A' we calculate the ni times multiplication {A' ) ni of the 
matrix A' . The non-zero elements of the matrix M ni ^ are 
equivalent to the non-zero elements of the matrix (A' ) ni . 
After determining the non-null upper off-diagonal blocks 
B^ k \ j < k of the matrix A^ satisfying the condition 
( |48| ) by such a process, the diagonal block B^' and the 
lower off-diagonal bocks B"*' , j > k of the matrix A^ 
is determined so that the matrix A<- n '' is symmetric and 



satisfies the condition (pf9|). ft should be noted that we 
could use a similar method to obtain the matrix A^ 
including the effect of long range interactions in the one 
dimensional model in the previous section. 

We restrict our consideration in the case that the 2x2 
matrix B«*' , which is the block element of the matrix A 
in the two dimensional system, is diagonalized. 



B 



m 
12 



B 



W 
21 







(59) 



This means that two components of the position of each 
particle do not interact with each other. 

Now we calculate the Lyapunov spectrum for such a 
matrix 

iM by using the formula (|34|). Figure ^| is the 
Lyapunov spectrum normalized by the maximum Lya- 
punov exponent in such a triangle lattice system in the 
case of Ni = 13, N 2 = 15 w (2/V3)Ni = 15.01110- •• 
and ni = 3. Here we chose the non-zero elements in the 
upper triangle of the matrix M ni ^ randomly from the pos- 
itive region (0.2tD, 1.8(D) with a (non-zero) real constant 
Co, and took the arithmetic average over this randomness. 
As in the one dimensional model including the long-range 
interactions, the Lyapunov spectrum is separated into the 
part changing smoothly and the part having a stepwise 
structure, which appears in a region of small positive 
Lyapunov exponents. 

A remarkable point in the two dimensional triangle 
lattice model, which the one dimensional models do not 
have, is the wide steps of the Lyapunov spectrum, espe- 
cially the step consisting of 4 points and the step consist- 
ing of 8 points shown in Fig. ^|. To give an explanation 
for these wide steps, we should notice that in this model 
we set-up a periodic boundary condition for each of the 
three directions of the triangle lattice. This boundary 
condition for each of the three directions can cause de- 
generacies in the Lyapunov spectrum like in the one di- 
mensional cases. Besides, in this model each site of the 
triangle lattice have two independent degrees of freedom 
corresponding to the two dimensionality of the particles, 
which can also cause degeneracy in the Lyapunov spec- 
trum due to the condition (|59|). Therefore we can get 
the step consisting of 4 points corresponding to each di- 
rection of the triangle lattice. Moreover two of the three 
directions (namely the two directions other than the hor- 
izontal direction in Fig. |^) have the same numbers of the 
sites in themselves, so they can cause a degeneracy in the 
Lyapunov spectrum. After all we get the step consisting 
of 4 points corresponding to the horizontal direction, and 
the step consisting of 8 points corresponding to the other 
two directions. This explanation is partly justified by the 
fact that the Lyapunov spectrum in the case of N\ = N 2 
has the step consisting of 12 points in the Lyapunov spec- 
trum, shown in the inset in Fig. ^, where we gave the 
averaged Lyapunov spectrum normalized by the maxi- 
mum Lyapunov exponent in the case of N\ = N 2 = 14 
and ni = 3. Here, except for the numbers N± and N 2 we 
used the same boundary condition and the same random- 
ness for non-zero elements of the matrix M ni ^ with in the 
case of N\ — 13 and N 2 — 15 in Fig. |(| This Lyapunov 
spectrum is shifted by —2 in the j-direction, so that the 
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positions of the zero Lyapunov exponents coincide with 
the cases of Ni = 13 and N2 — 15. One may also no- 
tice that the parts changing smoothly in the Lyapunov 
spectrum are almost indistinguishable in both the cases 
in the inset of Fig. 0. 



N,=13, N,=15 




Ni-13, N2=15 
Ni-14, N2=14 



FIG. 6. Lyapunov spectrum of the triangle lattice system 
in the case of Ni = 13, N 2 = 15 « (2/^)^1 and n ; = 3 
(the triangle dots). The inset: Enlarged graphs of the part 
including the stepwise structures in the Lyapunov spectra in 
the case of Ni = 13 and N2 = 15 (the triangle dots) and the 
case of Ni — N2 — 14 (the circle dots) for a comparison. The 
Lyapunov spectrum in the case of Ni = N2 = 14 is shifted 
by —2 in the j-direction, so that the positions of the zero 
Lyapunov exponents coincide in both the cases. 



Like in the one dimensional model, longer range inter- 
actions lead to a shorter region of stepwise structure in 
the Lyapunov spectrum in the two dimensional triangle 
lattice model. On the other hand, a numerical simula- 
tion of the many disk system showed that the region of 
the stepwise structure depends on the aspect ratio of the 
rectangular system This fact suggests that the long 
range interactions in the random matrix approach should 
depend on the aspect ratio of the rectangular system. 

It should be emphasized that such wide step consist- 
ing of 4 points or 8 points are actually found in the two 
dimensional deterministic chaotic system consisting of 
many hard-core disks numerically [B1J. Beside, we should 
also notice that in Fig. ^ (and Fig/g) a warp of the Lya- 
punov spectrum appears in a region of large Lyapunov 
exponents, which is also a characteristic of the system 
consisting of many hard-core disks. However there are 
some points which the triangle lattice model in this sec- 
tion cannot give enough explanation when compared with 
the numerically observed features of the Lyapunov spec- 
trum for the deterministic two dimensional hard disk sys- 
tem. First, in the numerical simulation of the hard-disk 
system each particle can interact with almost any other 
particle in the long time limit. This should correspond 
to a big number of ni w N1N2 in the triangle lattice 
model, but if we adopt such large number for ni then 
the stepwise structure disappears in the triangle lattice 
model. In this sense, the model in this section may cor- 
respond to the high density case, in which each particle 



mainly interacts with only a few particles surrounding it. 
Second, in the hard-disk system with a square shape the 
stepwise structure of the Lyapunov spectrum seems to be 
a repetition of the step consisting of 4 points and the step 
consisting of 8 points. On the other hand, in the triangle 
lattice model such a repetition of steps cannot be guar- 
anteed. As the third point, in the triangle lattice model 
in this section we adopt the boundary condition to make 
each of the three directions of the triangle lattice peri- 
odically, rather than the periodic boundary condition to 
make the up-side (the left-side) and the down-side (the 
right-side) of the square equivalently, which was adopted 
in the numerical simulation of the hard-disk system. 



VII. CONCLUSION AND REMARKS 

In this paper we have discussed the Lyapunov spec- 
trum in many particle systems described by a random 
matrix dynamics. We started from the many particle 
Hamiltonian mechanics without a magnetic field, and in- 
troduced Gaussian white random interactions between 
the particles. In such a system the dynamics of the 
tangent space is expressed by a Fokker-Planck equation, 
which leads to a direct connection between the positive 
(and zero) Lyapunov exponents and the time correlation 
of the matrix specifying the particles interactions. Us- 
ing this formula, we calculated concretely the Lyapunov 
spectra in one and two dimensional models satisfying the 
total momentum conservation with periodic boundary 
conditions. These models show a stepwise structure of 
the Lyapunov spectrum in the region of small positive 
Lyapunov exponents, which is robust to a perturbation in 
matrix elements of the particle interaction matrix. The 
long range interactions between the particles lead to a 
clear separation into a part exhibiting stepwise structure 
and a part changing smoothly. The part of the Lyapunov 
spectrum containing stepwise structure is clearly distin- 
guished by a wave like structure in the eigenstates of the 
particle interaction matrix. In the two dimensional model 
we got wider steps in the Lyapunov spectrum than in the 
one dimensional models, especially the steps consisting of 
4 points and the steps consisting of 8 points. These wide 
steps in the Lyapunov spectrum have already been shown 
numerically in a deterministic chaotic system consisting 
of many hard-core disks. 

One of the important simplifications in this random 
matrix approach is that in this approach we do not have 
to refer to the phase space dynamics in order to deter- 
mine the tangent space dynamics any more. In general, 
the matrix R(t) appearing in Eq. (0) can depend on the 
phase space dynamics, so this kindof separation of the 
phase space dynamics and the tangent space dynamics 
does not be allowed in deterministic chaotic systems. 

As emphasized in the discussion of the two dimensional 
model, a lattice picture is useful to make a model in the 
random matrix approach. Concerning this point it may 
be interesting to note that some works for the Lyapunov 
spectra for many particle systems indicated a similarity 
between solid state phenomena and the behavior of the 
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Lyapunov spectra for many particle systems in a fluid 
phase 1 16 37]. 

We have regarded the random matrix dynamics in this 
paper as an imitation of the deterministic chaotic dy- 
namics, and reproduced some characteristics of chaotic 
systems, especially the stepwise structure of the Lya- 
punov spectra. However we must not forget that there 
are some differences between the chaotic dynamics and 
the random matrix dynamics in this paper. For example, 
in the random matrix dynamics the movement of parti- 
cles are not deterministic but stochastic, so that the zero 
Lyapunov exponents arising from the initial infinitesimal 
perturbation along the orbit in the deterministic chaos do 
not appear in the random matrix dynamics. (Note that 
we got only the zero Lyapunov exponents correspond- 
ing to the total momentum conservation in the models 
discussed in this paper.) We should also mention that 
in the random matrix approach there is an additional 
statistical average over the randomness of the interac- 
tions of particles, which does not exist in the determin- 
istic dynamics. This causes vagueness in the definition 
of the Lyapunov exponents. In this paper we introduced 
the Lyapunov exponent as the time average of the log- 
arithm of the randomness-averaged time evolution of a 
neighboring trajectory. This definition allows us to get a 
simple connection between the Lyapunov exponents and 
the time correlation of the interaction matrix as shown 
in Eq. (26). However this definition of the Lyapunov 



exponent is not proper to discuss the negative Lyapunov 
exponents and hence the pairing rule for the Lyapunov 
spectrum. On the other hand, we could adopt the defi- 
nition of the Lyapunov exponent as the time average of 
the randomness-average of the logarithm of time evolu- 
tion of a neighboring trajectory. This definition should 
be proper to discuss the negative Lyapunov exponents. 
The comparison of these two definitions of the Lyapunov 
exponent is an unsettled problem. 

One may regard the random matrix approach using 
the master equation in this paper as one of the few ana- 
lytical approaches to calculate the Lyapunov exponents. 
The other statistical and analytical approach for the Lya- 
punov exponents is the kinetic theoretical approach of 
Refs. p8|j3l| . In this approach the positive Lyapunov ex- 
ponents are calculated using a Lorentz-Boltzmann equa- 
tion, while the negative Lyapunov exponents are calcu- 
lated using an " anti-Lorentz-Boltzmann equation" where 
the collision operator has the opposite sign to the ordi- 
nary Lorentz-Boltzmann equation. However, so far this 
kinetic theoretical approach can only provide the maxi- 
mum Lyapunov exponent and the Kolmogorov-Sinai en- 
tropy for dilute gases. 

To improve the random matrix approach used in this 
paper it is essential to know the statistical information for 
the interaction matrix R(t). It has already been shown 
numerically that the time average of the matrix R(t) is 
almost null in the system consisting of many hard-core 
disks EfJ. This result justifies the condition (^) in the 
case ofn = 1. A similar investigation of the correlation 
of the matrix R(t) in deterministic many particle chaotic 
systems is one of the important future problems. 
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APPENDIX A: MASTER EQUATION FOR THE 
TANGENT SPACE 



In this appendix we derive the Fokker-Plank equation 
(O) for the tangent vector space. Using the Kramers- 
Moyal expansion the dynamics of the probability density 
p(ST, t) is given by 



dp(ST, t) 
di 



oo 2N 2N 2N 

£££•••£(-!)" 

n=ljl=lj2 = l Jn=l 

d5r h dsr h ---dsr jn 



where ^j™j 2 ...j„(^I\ t) is defined by 

x[8r h (t + s )-6r h (t)} 
■■■[sr jn (t + s)-sr jn (t)} 



(A.l) 



sr(t)=sr 



(A.2) 



and ST-jU) is the j-th component of the tangent vector 

5v(t) m. 

It follows from Eqs. @, (§), (|) and (§) that 



ST(t + s) - ST(t) 

r t+s 

J i dr L(t) 



£ 



T exp 

t+s 



<iT n J <lr„ : 



- 1 ST(t) 



dr\ 



xJL{T n )JL(T n ^) ■ ■ ■ JL{ Tl )5T(t) 

and 

JL(T. n )JL(T n ^ 1 )---JL{T 1 ) 



(A.3) 



I/m 
R(ti) 



(i) 
i 



$ 





(2) 

i 



for n = 1 



for n = 21,1 = 1,2,- 



(1) 1 ' for n = 2/ + l,Z = 1,2,- 



(+1 



(A.4) 
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where $ ; , j — 1,2 are defined by 

$W = i?(T 2; _ 1 )i?(r 2; _ 3 ) • • • R(Ti)/m l (A.5) 

$ ; (2) = i?(r 2i )i?(r 2; _ 2 ) ■ • • R(T 2 )/m l . (A.6) 
By using Eqs. (g), @, ( |A~2| ), ( |A~3| ) and (fO|) we obtain 

s (1) (5r,i) 

= l^([ST(t + s)-8T(t)})\ sm=ST 
1 f t+s 

= lim - / dr (JL(t)) 5T 
s^O s J t 

= (5^, <^ 2 ■ • • , Sp N , 0,0, • • • , 0) T /m (A.7) 

H( 2 )(«5r,t) = (Hg)(,5r,f)) 



= iim-([ar(t + a )-ar(f)] 

x [ST(t + s)-ST{t)] J 

1 f t+s f t+s 
= lim — / d,K / dr 
2s J t J t 



JL(k)5T5T t [JL{t) 





*( 2 H<5q) 

where * (2) ((5q) = (*^(5q)) is defined by 

Q=l (9=1 



«r(t)=<5r 



(A. 



(A.9) 



Here the only non zero contributions come from the n = 1 
term of Eq. (A. 3). For general n, the number of delta 
functions from Eq. (10) must be only one less than the 
number of time integrals, to give a non zero contribution. 
It is straightforward to show that this never happens for 
n > 1. The terms including S^ ...^ [ST, t), n = 3, 4, • • ■ 
in the right-hand side of Eq. ( A.l) are negligible because 
of the Gaussian white properties (^) and ([l0]) of the ran- 
dom matrix i?(£). Using this fact and Eqs. (|l2|), (A.l), 
(A.7) and (A. 8) we obtain the Fokker-Planckequation 



APPENDIX B: LYAPUNOV EXPONENTS IN 
THE RANDOM MATRIX APPROACH 

In this appendix we derive the expression ( B 6 [) for the 
Lyapunov exponents from the definition ([251)7 Trans- 
forming Eqs. (|~~), ( |l~ ) and ( ~(i| ) into the equations for 

the quantities T^2(t) we obtain 



dt 



7TJ ■* 



(B.l) 
(B.2) 



dt 



N N N N 

a=l 0=1 M=l "=1 
AT iV W iV 

a=l (5=1 ii=l v-1 

N N N N 

x E E EEw^^w 



Q!' = l /3' = 1 /i' = l v' = l 



(B.3) 



where we used 



(B.l), (B.2) and 



EqMO), ©, (||) and Equations 
d (pT3neadto 



dt 3 

'33 
dt 2 



,2 « 



dT ^ = lf(3) (t) 

m 2 JJ 

dT 5?w_ 2- (2) 



77/ 



(B.4) 
(B.5) 
(B.6) 



It is noted that Eq. (B.4) is the differential equation only 



for the quantity T^(t) 

If the quantity Ajj is zero, then the quantity (t) is 
a bilinear function of time t so that the Lyapunov expo- 
nent defined by Eq. (^5|) gives zero, namely Eq. ( pq) is 
correct in this case. In the case of Ajj ^ 0, notingthat 
the functions exp[(2A w /m 2 ) 1 / 3 t-exp(27rA:i/3)], k = 0, 1, 2 
are special solutions of Eq. (B.4), we obtain the general 
solution 



-Re | (nf ) - inf ] j eq> 



/Cji • exp — 7ri 



f^ iJ exp (/Cj-t) 



fit. 2 ) cos h^-Kjt 



j 



flf sin(f^ 



x exp ( ——K-jt 



(B.7) 



of Eq. (B.4) for the real function Y^?(t) where /Cj is 
defined by 



(2A 



jj 



\ m 2 



1/3 



(B- 
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Here £lf\ k = 1,2,3 are real constants and are connected 
to the initial condition as 

2T^(0) 







Tg>(0) + 



2T^(0) 



mK,j (m/Cj) 5 



n?> = 4 



n< 3) = 



tg(Q) 

{mJCj) 2 



llCjy/3 



tg } (0) 
m/C, 



(B.9) 
(B.10) 
(B.ll) 



by u sing Eqs. 
Eq. (fjj) into E"q~ 



B.6 ) and (B.7). By substituting 
we obtain 



x - ^ 



A; 



(2m) 2 



1/3 



(B.12) 



+0jvj0j2 exp(-27ri(j - l)k/N) 
+6(N-i)j0n cxp(-27ri(j + l)k/N) 
+O N] 0]N cxp(-27ri(j - TV + l)fc/JV) 
+0ij0ji cxp(-2vri(j + N — l)k/N)] 



where 0jfc is defined by 



1 in j > k 
in j < k 



(C.2) 



(C.3) 



and satisfies the relations 9jk&kj = <5jfc and 0/v„0 n ; + 
0(/_i) n 0ni = 1 for any integer {2, 3, • • • , TV] and any 

integer nS {1, 2, • • • , AT}. Here the A-th eigenvalue a^ 1 of 
the matrix is given by 



namely Eq. ( 1261) . 

Equations ( p.5[) and ( B.6 ) imply that the quantities 



(o) 



-2w 1 



2?rfc 



(i) and (i) are given by taking the first and sec- 
ond derivative of the quantity (t) with respect to the By using the formula ( |34[ ) for the eigenvalues a 
time t, respectively. This fact leads to the relation 

H1 \t) l Tg>(t) 

hm — In 

t— >oo 2i 



(C.4) 



,(o) 



1 T 
lim — ln-ryf- 

t^oo2t T i)( ) 



T 

n 

, y(3) 



lim — In ■ 

t^oo 2i 



(2) (0) 



(B.13) 



Tg>(0) 



et the same Lyapunov exponents as 
!9) through the equations lim t ^ +co 

(2/)- 1 HT^(t))/f5; ) (0)],A 



Namely we 
in Eq. 



2,3. 



APPENDIX C: LYAPUNOV EXPONENTS FOR 
ONE DIMENSIONAL MODELS WITH THE 
NEAREST NEIGHBOR INTERACTIONS 



In this appendix we give the derivation of Eqs. ( pq ) 
and (||) from Eqs. (g) and @. We also show Eqs. 
© and (H) in the case of A = yl (0) + e^ 1 ' . 

We define the discretized Fourier-transform matrix 
F = (F 3k ) by 



1 



exp 



-2ni ■ 



A:/ 



given by Eq. (0.4) we obtain Eqs. (|38[). The eigenval- 
ues of the matrix A^ have degeneracies because of the 
relation a^_j = a^ ' in j < N/2, so we obtain Eq. (|39|). 
Equation (C.2) implies that the matrix is diago- 



nalized by using the matrix F: (F^ A^ F)jk 



,(0); 



(B.14) rp ne eigenvectors xf^ of the matrix A^ corresponding 



to the eigenvalue is represented as 



xf = (F lj ,F 2j ,...,F Nj f 



(C5) 



so that we have the relation A^x^ = a^XnK This 
set of eigenvectors satisfies the completeness condition 

and is normalized, namely 5Z^ =1 xi°^Xa t = I anc ^ 
x (o) t (o) = , 

Second, we expand the j-th eigenvalue aj of the ma- 
trix A = A^ + eA^ in the small parameter e, namely 
aj = a^ "* + ea~- + • • •, and consider the first order correc- 
tion £dj by using the well-known perturbation theory 
for a degenerate system. We introduce the eigenstate Xj 
of the matrix A corresponding to the eigenvalue Oj , and 

(oh 



(C.l) expands it with the complete set {x\ '}i of the vectors 



which is an unitary matrix: F*F — FF^ — I with the 
superscript f representing the Hermitian conjugate of the 
matrix. 

First, it follows from Eqs. @ and (|oT|) that 

N 

= -^[-2exp(-27rijfe/^) 



JV 

E< 

a=l 



(C6) 



with the constant Cj a = x& ^Xj. The relation Ax,,- 



OjXj is translated into the equation 



/v 



[(«. 



-aj)5 ka +eW ko 



(C.7) 
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for the coefficients {cjk}j,k and the eigenvalue a,j with in the case of j < N/2, where we used the relation 
the quantity 



\W. 



W ik = ^A^h 



.(0) 



We expand the coefficient Cj k by the small parameter e: 
•. Now we calculate the first or- 



„(0) 



Cjk = C'jk + £ C jk 

der corrections ectj 1 ' and ea^_^ of the eigenvalues aj and 
djv-j; respectively, which have a degeneracy in the 0-th 
order of the parameter e. For this purp ose, instead of the 
coefficients {cjk}j,k appearing in Eq. ( p.7p it is enough 

to consider the 0-th order coefficients {cf k }j. k , which are 

zero except for c 



(i) 



(o) 

33 ' 



This leads to the equation 



(o) (o) 

C 3{N-3V C {N-j)j 



(0) 



Det 



-a« + Wi- 



tt', 



(N-j)j 



-aP> 



w W-j) 

+ W {N _ j)(N _ 3) 



= (C.9) 



3(N~3)\ 



XjCij /u}\ and numbered so that we obtain 



(C.8) (a) 1 ' - a$_jWoy j > 0. By using Eq. (|C.12| ) we also 



,(o) 



obtain ffl^JL = Xocjyyj/' 4 ' ' n ^ ne case °^ ^ nc num ber N to 



,(0) 



be even, and 



l N 



for the first order correction to the 



non-degenerate eigenvalues of the matrix . By apply- 
ing the formula @) to the case of a,j — '+eop +0(e 2 ) 



with the quantity a,- given by Eq. (C.14) or (C.15), we 



obtain Eqs. (^lj) and d42j) 



APPENDIX D: LYAPUNOV EXPONENTS FOR 
ONE DIMENSIONAL MODELS WITH LONG 
RANGE INTERACTIONS OF A FIXED 
STRENGTH 



for whose solutions give the eigenvalues a*, and 

a (1) 

a N-j 



We can solve Eq. flC.9|) and obtain 



W (N~ 3 )(N-j) 




W {N ^ j){N -j) 



+ \ W 3(N-3)\ , 



(CIO) 

noting the relation VFfcj = with the superscript * 
representing the complex conjugate of the complex num- 
ber. The eigenvalue a$ has no degeneracy so we sim- 
ply have a $ — Wnn- If the number N is even, then 
the eigenvalue a^/ 2 has no degeneracy, so we also have 



"N/2 ~ w ( N /2)(N/2) in this case. 
By using Eqs. ©, (|(Tl] ), flcj ) and we obtain 

= - [1 - exp {-2mk/N)] 
x [l-exp(27r^/7V)]iV- 1 



x > Xa exp [27ria(fc - l)/N] . 



a=l 



(C.ll) 



In this appendix we calculate the Lyapunov spectra 
for the systems with long range interactions of a fixed 
strength. 

First we consider the case described by the matrix 
A = It"'] = defined by 



+8j(k+N-t) + 3k(j+N~l)] 



(D.l) 



with ni < N/2 and a (non-zero) real constant uj. It 
should be noted that the matrix A\ ni ' given by Eq. 
(Eft) is attributed into this matrix A^ in the case of 

a j = Q. The matrix A^ is diagonalized by the ma- 



trix F defined by Eq. (C.l), and we obtain the rela 



tion (F^A^F) jk = a [ k l] 8 jk with the eigenvalue af' J = 

— 2tD [ni — Yl^Li cos {^jl/N)] . By applying the formula 

( p4| ) to the eigenvalue a,j = 
exponent 



we obtain the Lyapunov 



Especially we derive 



W kk = - |1 - exp (-2mk/N)\ 2 X o 

(°) w 
= — a k = W(JV-fc)(AT-fc) 



k(N-k) 



= - [1 - exp (-2™fc/iV)] 2 Xfe 



(C.12) 



(C.13) 



f rom E qs. 
(|012|) and 



( C.4| ) and (|C.llj) . By substituting Eqs. 
lMlnto Eq. (p.l0|) we obtain 



(0) 

— IXo 

l N-j 



\Xj\ 



(xo - IXil 



(C.14) 
(C.15) 



2/3 



/ 2 71^7 

n; — > cos 



i=i 



2/3 



(D.2) 



The Lyapunov exponents given by Eq. ( |D.2| ) satisfies 
the relation Xj = Ajv_j in j < N/2, so the Lyapunov 
spectrum of this system has a stepwise structure. 

Second we consider the case that each particle inter- 
acts with all the other particles with the same strength. 

This is described by the matrix A = A^^ = {Af k /2] ), 
which is defined by 



.4 



[N/2] 
jk 



2w(l - NS jk ), 



(D.3) 
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namely the matrix whose off-diagonal elements are n on- 
zero and equal with each other. The matrix (D.3) is 
diagonalized as 



^A [N/2] F) jk = -2NQ(1 - 8 jN )S. 



jk 



(D.4) 



by using the matrix F defined by Eq. (CJ), so the eigen- 
values of the matrix A^ N / 2 ^ are —2Nu> and 0. By substi- 
tuting the eigenvalues of the matrix A — A^- N ^ 2 ^ in the 
formula (34) we obtain the Lyapunov spectrum as 



A b1 = 



NlD i 2 / 3 
m I 

o 



in j 



1,2,- 

N, 



■N - 1 



(D.5) 



namely the shape of positive Lyapunov spectrum in this 
system is just in a straight horizontal line. It should be 
noted that the quantity Q may depend on the number N 
of the particles in general. Therefore the consideration in 
this appendix is not enough to discuss the particle num- 
ber dependence of the maximum Lyapunov exponent. 
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